In  this  report  we  investigate  the  solution  of  boundary  value  problems  on  polygonal 
domains  for  elliptic  partial  differential  equations.  We  observe  that  when  the  problems 
are  formulated  as  the  boundary  integral  equations  of  classical  potential  theory,  the 
solutions  are  representable  by  series  of  elementary  functions.  In  addition  to  being 
analytically  perspicuous,  the  resulting  expressions  lend  themselves  to  the  construction 
of  accurate  and  efficient  numerical  algorithms.  The  results  are  illustrated  by  a  number 
of  numerical  examples. 
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1  Introduction 


In  classical  potential  theory,  elliptic  partial  differential  equations  (PDEs)  are  reduced  to 
integral  equations  by  representing  the  solutions  as  single-layer  or  double-layer  potentials 
on  the  boundaries  of  the  regions.  The  densities  of  these  potentials  satisfy  Fredholm 
integral  equations  of  the  second  kind. 

There  are  three  essentially  separate  regimes  in  which  such  boundary  integral  equa¬ 
tions  have  been  studied.  In  the  first  regime,  the  boundary  of  the  region  is  approximated 
by  a  smooth  curve.  It  is  known  that  if  the  curve  is  smooth,  then  the  kernel  of  the  in¬ 
tegral  equation  is  smooth  as  well  (see,  for  example,  [21]).  The  existence  and  uniqueness 
of  the  solution  follows  from  Fredholm’s  theory,  and  the  integral  equations  can  be  solved 
numerically  using  standard  tools  (see,  for  example,  [11]). 

In  the  second  regime,  the  boundary  of  the  region  is  approximated  by  a  curve  with 
perfectly  sharp  corners.  In  this  regime,  the  kernel  of  the  integral  equation  has  singulari¬ 
ties  at  the  corners,  and  the  existence  and  uniqueness  of  the  solution  in  the  L2-sense  is  also 
known  (see,  for  example,  [23]).  The  behavior  in  the  vicinity  of  the  corners  of  the  solu¬ 
tions  of  both  the  integral  equations  and  of  the  underlying  differential  equation  have  been 
the  subject  of  much  study  (see  [25],  [15]  for  representative  examples).  Comprehensive 
reviews  of  the  literature  can  be  found  in  (for  example)  [19],  [10]. 

In  the  third  regime,  the  assumptions  on  the  boundary  are  of  an  altogether  different 
nature.  It  might  be  a  Lipschitz  or  Holder  continuous  curve,  or  a  fractal,  etc.  While 
during  the  last  fifty  years,  such  environments  have  been  studied  in  great  detail  (see,  for 
example,  [13],  [23],  [4],  [6],  [5],  [14],  etc.),  they  are  outside  the  scope  of  this  paper. 

This  paper  deals  with  the  very  special  case  of  polygonal  boundaries,  and  is  based 
on  several  specific  analytical  observations.  These  observations  are  summarized  in  the 
following  two  subsections,  and  are  discussed  in  detail  in  sections  4  and  5. 
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1.1  The  Neumann  Case 

Suppose  that  7:  [—1, 1]  — >  R2  is  a  wedge  in  R2  with  a  corner  at  7(0),  and  with  interior 
angle  not.  Suppose  further  that  7  is  parametrized  by  arc  length,  and  let  v{t)  denote 
the  inward-facing  unit  normal  to  the  curve  7  at  t.  Let  Y  denote  the  set  7Q— 1, 1]).  By 
extending  the  sides  of  the  wedge  to  infinity,  we  divide  R2  into  two  open  sets  JR  and  O2 
(see  Figure  1). 

Let  4>:  R2  \  T  — *  R  denote  the  potential  induced  by  a  charge  distribution  on  7  with 
density  p:  [—1, 1]  — >■  R.  In  other  words,  (f)  is  defined  by  the  formula 

<i>(x)  =  J  log(||7(i)  —  x\\)p(t)  dt,  (1) 

for  all  x  G  R2  \  T,  where  ||  •  ||  denotes  Euclidean  distance.  Suppose  that  n  is  a  positive 
integer,  and  let  |~ x]  denote  the  smallest  integer  n  such  that  n  >  x,  and  |_r|  denote  the 
largest  integer  n  such  that  n  <  x,  for  all  real  x.  Suppose  further  that  p  is  defined  by  the 
formula 


\n/2\  \n/2\  \n/2\  \n/2\ 

1  2 i — 1  -1  x  ^  ,  /  \  1  1  2i  1  x  x  /  \  1  1  2i  1  1  x  ^  1  1  2i 


p(t)  =  ^  hi-i\t\  “  1  +  Y  &2*sgn(f)|f|«  1  +  Y  c2i-i  sgn(f)|f|  2'-a  1  +  Y  c2 i\t\2~a  \ 
i=l  1=1  1=1  1=1 

(2) 

for  all  —  1  <  t  <  1,  where  b\,  b-2,  ■  ■  ■ ,  bn  and  ci,  c2, . . . ,  cn  are  real  numbers,  and 


—  1  if  x  <  0, 
sgn(x)  =  <(  0  if  x  —  0, 

1  if  x  >  0, 

for  all  real  x.  Suppose  finally  that  g :  [—1, 1]  — *  R  is  defined  by  the  formula 

,  .  ..  d (j>(x) 

9(t)=  hrn  — — 

x—>'y(t)  OlSyt) 


(3) 


(4) 


for  all  —  1  <  t  <  1,  i.e.  g  is  the  limit  of  the  normal  derivative  of  integral  (1)  when 
x  approaches  the  point  7 (t)  from  outside.  Then  g  is  smooth  on  each  of  the  intervals 
[—a,  0],  [0,  a]  for  any  0  <  a  <  1  (see  Figure  1).  In  other  words,  given  a  charge  distribution 
p  of  the  form  (2),  the  normal  derivative  of  the  resulting  potential  (1)  on  Y  (in  the 
sense  (4))  is  smooth. 

Conversely,  for  any  sufficiently  smooth  g,  there  exists  a  charge  distribution  p  of  the 
form  (2)  such  that  the  normal  derivative  of  the  resulting  potential  (1)  on  Y  is  equal  to 
g,  to  high  precision. 
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1.2  The  Dirichlet  Case 


Suppose  that  7:  [—1, 1]  — >  R2  is  a  wedge  in  R2  with  a  corner  at  7(0),  and  with  interior 
angle  not.  Suppose  further  that  7  is  parametrized  by  arc  length,  and  let  v{t)  denote 
the  inward-facing  unit  normal  to  the  curve  7  at  t.  Let  Y  denote  the  set  7Q— 1, 1]).  By 
extending  the  sides  of  the  wedge  to  infinity,  we  divide  R2  into  two  open  sets  IR  and  122 
(see  Figure  1). 

Let  (j):  R2  \  T  — *  R  denote  the  potential  induced  by  a  dipole  distribution  on  7  with 
density  p\  [—1,1]  — *  R.  In  other  words,  c j>  is  defined  by  the  formula 


Hx) 


(z/(f),7(f)  -x) 
||7(f)  -x\\2 


pit )  dt, 


(5) 


for  all  x  G  R2  \  T,  where  (•,  •)  denotes  the  inner  product  and  ||  •  ||  denotes  Euclidean 
distance.  Suppose  that  n  is  a  positive  integer,  and  let  |"x]  denote  the  smallest  integer  n 
such  that  n  >  x,  and  [rrj  denote  the  largest  integer  n  such  that  n  <  x,  for  all  real  x. 
Suppose  further  that  p  is  dehned  by  the  formula 


\n/2\  \n/2\  \n/2]  _  \n/2\ 

pit)  =  Y  h2i-l\tf^  +  Y  ^2jSgn(f)|f|«  +  Y  C2i-lSgn(f)|f|^  +  Y  C2i\t\^- 
2—1  2—0  2=1  2=0 


(6) 


for  all  —  1  <  t  <  1,  where  b0,bi, . . .  ,cn  and  Co,  Ci, . . . ,  cn  are  real  numbers,  and  sgn(a;)  is 
dehned  by  (3).  Suppose  finally  that  g :  [—1, 1]  — *  R  is  dehned  by  the  formula 

9it)  =  lim  Mx),  (7) 

for  all  —1  <  t  <  1,  i.e.  g  is  the  limit  of  integral  (5)  when  x  approaches  the  point  7 (f) 
from  inside.  Then  g  is  smooth  on  each  of  the  intervals  [— a,  0],  [0,a]  for  any  0  <  a  <  1 
(see  Figure  1).  In  other  words,  given  a  dipole  distribution  p  of  the  form  (6),  the  resulting 
potential  (5)  on  Y  (in  the  sense  (7))  is  smooth. 

Conversely,  for  any  sufficiently  smooth  g,  there  exists  a  dipole  distribution  p  of  the 
form  (6)  such  that  the  resulting  potential  (5)  on  Y  is  equal  to  g,  to  high  precision. 


1.3  The  Procedure 

Recently,  progress  has  been  made  in  solving  the  boundary  integral  equations  of  potential 
theory  numerically  (see,  for  example,  [12],  [3]).  Most  such  schemes  use  nested  quadratures 
to  resolve  the  corner  singularities.  However,  the  explicit  representations  (2),  (6)  lead  to 
alternative  numerical  algorithms  for  the  solution  of  the  integral  equations  of  potential 
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theory.  More  specifically,  we  use  these  representations  to  construct  purpose-made  dis¬ 
cretizations  which  accurately  represent  the  associated  boundary  integral  equations  (see, 
for  example,  [18],  [16],  [24]).  Once  these  discretizations  are  available,  the  equations  can 
be  solved  using  the  Nystrbm  method  combined  with  standard  tools.  We  observe  that 
the  condition  numbers  of  the  resulting  discretized  linear  systems  closely  approximate  the 
condition  numbers  of  the  underlying  physical  problems. 

Observation  1.1  While  the  analysis  in  this  paper  applies  only  to  polygonal  domains, 
a  similar  analysis  carries  over  to  curved  domains  with  corners.  A  paper  containing  the 
analysis,  as  well  as  the  corresponding  numerical  algorithms  and  numerical  examples,  is 
in  preparation. 

Observation  1.2  In  the  examples  in  this  paper,  the  discretized  boundary  integral  equa¬ 
tions  are  solved  in  a  straightforward  way  using  standard  tools.  However,  if  needed,  such 
equations  can  be  solved  much  more  rapidly  using,  for  example,  [9]. 

The  structure  of  the  paper  is  as  follows.  In  Section  2,  we  introduce  the  necessary 
mathematical  preliminaries.  Section  3  contains  the  primary  analytical  tools  of  the  pa¬ 
per.  In  sections  4  and  5,  we  investigate  the  Neumann  and  Dirichlet  cases  respectively. 
In  Section  6,  we  briefly  describe  a  numerical  algorithm  and  provide  several  numerical 
examples. 

2  Mathematical  Preliminaries 

2.1  Boundary  Value  Problems 


7(ss) 


Figure  2:  A  curve  in  R2 

Suppose  that  7  :  [0,  L]  — y  R2  is  a  simple  closed  curve  of  length  L  with  n  corners  at  the 
points  0  =  si  <  S2  <  ■  ■  ■  <  sn  <  sn+i  =  L  .  Suppose  further  that  7  is  parameterized  by 
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arc  length  and  oriented  counterclockwise.  We  denote  the  interior  of  7  by  and  denote 
the  boundary  of  by  T.  Suppose  that  7  is  analytic  on  the  intervals  (sj,  si+i)  for  each 
i  =  1,2, ...  ,n,  and  that  the  interior  angle  at  any  corner  of  7  is  between  0  and  2n.  We 
denote  the  normalized  internal  normal  to  7  at  i  £  [0,  L\  by  v(t).  Letting  g  be  a  function 
[0,  L\  — y  R,  we  will  consider  the  following  problems. 

Exterior  Neumann  problem:  find  a  function  0:  R2  \  — >  R  such  that 

V20(x)  =  0  for  x  G  R2  \  H,  (8) 

=  g(t )  for  t  G  [0,  L\.  (9) 

x— py(t)  OlSyt) 
xElR?\fl 

Interior  Dirichlet  problem:  find  a  function  0:  — »  R  such  that 

V20(x)  =  0  for  a;  G  (10) 

lim  0(x)  =  g(t)  for  t  G  [0,  L].  (11) 

x— >7(t) 


For  every  (7  G  L2([0,L]),  the  exterior  Neumann  problem  and  the  interior  Dirichlet  prob¬ 
lem  have  unique  solutions  (see,  for  example,  [13]). 


2.2  Integral  Equations  of  Potential  Theory 


In  classical  potential  theory,  boundary  value  problems  are  solved  by  representing  the 
function  0  by  integrals  of  potentials  over  the  boundary.  The  potential  of  a  unit  charge 
located  at  x0  G  R2  is  the  function  0°o :  R2  \  x0  — >  R,  defined  via  the  formula 

^x0(x)  =  Mlk-^oll),  (12) 

for  all  x  G  R2  \  x0,  where  ||  •  ||  denotes  Euclidean  distance.  The  potential  of  a  unit 
dipole  located  at  x$  G  R2  and  oriented  in  direction  h  G  R2,  \\h]\  =  1,  is  the  function 
h :  R2  \  Xq  — >  R,  defined  via  the  formula 


0 


1 

xo,h 


(x) 


{h,x 0  -  x) 

ll^o  -  a;||2 


(13) 


for  all  x  G  R2  \  x0,  where  (•,  •)  denotes  the  inner  product. 

The  potential  created  by  a  charge  distribution  with  density  p:  [0,  L\  — >  R  on  T  is 
given  by  the  formula 


00*0 


(14) 


for  any  x  G  R2  \  T,  and  is  called  a  single  layer  potential.  The  potential  created  by  a 
dipole  distribution  with  density  p :  [0,  L\  — >  R  on  T  is  given  by  the  formula 

(p(x)=[  ip^{t)At)(x)p(t)dt,  (15) 

Jo 

for  any  ieR2\r,  and  is  called  a  double  layer  potential. 

2.2.1  Reduction  of  Boundary  Value  Problems  to  Integral  Equations 

The  following  two  theorems  reduce  the  boundary  value  problems  of  Section  2.1  to  bound¬ 
ary  integral  equations.  They  are  found  in,  for  example,  [23]. 

Theorem  2.1  Suppose  that  p  G  L2([0,L]).  Suppose  further  that  g :  [0,  L]  — »  R  is  defined 
by  the  formula 

g(s)  =  -7T p(s)  +  [  ^7(a),„(i)(7 (t))p(t)dt,  (16) 

Jo 

for  any  s  G  [0,  L\.  Then  g  is  in  L2([0,L]),  and  the  solution  to  the  exterior  Neumann 
problem  with  right  hand  side  g  is  given  by  (If).  Moreover,  for  any  g  G  L2([0,L]), 
equation  (16)  has  a  unique  solution  p  G  L2([0,L]). 

Theorem  2.2  Suppose  that  p  G  L2([0,L]).  Suppose  further  that  g :  [0,  L\  — >  R  is  defined 
by  the  formula 

g(s)  =  -np(s)  +  f  ^(t)iI/(t)(7 (s))p(t)dt,  (17) 

Jo 

for  any  s  G  [0,  L].  The  g  is  in  L2([0,L]),  and  the  solution  to  the  interior  Dirichlet 
problem  with  right  hand  side  g  is  given  by  (15).  Moreover,  for  any  g  G  L2([0,L]), 
equation  (17)  has  a  unique  solution  p  G  L2([0,L]). 

Observation  2.1  Equation  (16)  is  the  adjoint  of  equation  (17). 

Observation  2.2  Suppose  that  the  curve  7:  [0,  L\  — >  R2  is  not  closed.  We  observe  that 
if  p  G  L2([0,L]),  and  g  is  defined  by  either  (16)  or  (17),  then  g  G  L2([0,L]).  Moreover, 
if  g  G  L2([0,L]),  then  both  equations  (16)  and  (17)  have  unique  solutions  p  G  L2([0,L]). 
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2.2.2  Properties  of  the  Kernels  of  Equations  (16)  and  (17) 

The  following  theorem  shows  that  if  a  curve  7  is  sufficiently  smooth,  then  the  kernels  of 
equations  (16)  and  (17)  are  bounded.  It  is  found  in,  for  example,  [20]. 

Theorem  2.3  Suppose  that  7:  [0,  L\  — >■  R2  is  a  curve  in  R2  that  is  parameterized  by 
arc  length,  and  that  is  not  necessarily  closed.  If  7  is  C 2  in  a  neighborhood  of  a  point  s, 
where  0  <  s  <  L,  then 

\  =  =  -\k(Ks)i  (18) 

where  k\  [0,  L]  — >  R  is  the  signed  curvature  of'y. 

The  following  corollary  shows  that  if  7  is  analytic,  then  the  kernels  of  equations  (16) 
and  (17)  are  also  analytic. 

Corollary  2.4  Suppose  that  7:  [0,  L\  — >■  R 2  is  a  curve  in  R2  that  is  parameterized  by 
arc  length,  and  that  is  not  necessarily  closed.  Ifj  is  analytic  in  a  neighborhood  of  a  point 
s,  where  0  <  s  <  L,  then 


^WM,)(7W),  (19) 

^7(t)Mt)(T'(s))>  (20) 

are  analytic  functions  of  t  in  a  neighborhood  of  s. 

Proof.  We  observe  that  there  is  some  neighborhood  V  of  s  such  that 


(21) 

^7(t)M*)(T'(s))>  (22) 


are  analytic  functions  of  t  on  K\{s}.  Since,  by  Theorem  2.3,  these  functions  are  bounded 
at  s,  it  follows  that  they  are  analytic  at  s.  ■ 


When  the  curve  7  is  a  wedge,  the  kernels  of  equations  (16)  and  (17)  have  a  particularly 
simple  form,  which  is  given  by  the  following  lemma. 


Lemma  2.5  Suppose  7:  [—1, 1]  — >  R2  is  defined  by  the  formula 


l(f) 


—t  ■  (cos(7tq;),  sin(7ra))  if  —  1  <  t  <  0, 
(7,0)  if0<t<l, 


(23) 
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shown  in  Figure  3.  Then ,  for  all  0  <  s  <  1, 

{isin(7ra)  _ 

— - - - - -  if  —  1  <  t  <  0, 

s2  +  t2  +  2stcos(na) 

0  if0<t<l, 


and,  for  all  —  1  <  s  <  0, 


[  s2  +  t2  +  2 st  cos(7ra) 
Proof.  Suppose  that  0  <  s  <  1  and  0  <  t  <  1.  Then, 


if  —  1  <  t  <  0, 
if0<t<l. 


i  „  (v(s),ry(s)  —  7 (t)) 

=  ||7(s)  _7(t)||2 

((0,l),(5,0)-(t,0)) 

Is  -  A2 

=  0. 

Now  suppose  that  0  <  s  <  1  and  —  1  <  t  <  0.  Then, 


,i  ( r\/( \  -  Ms),7W-7W) 
^(s)As)h(t))  ||7(s)  -7(t)||2 

((0, 1),  (s,  0)  +  t(cos(7ro;),  sin(7ra))) 
(s  +  t  cos(7TQ;))2  +  (tsin(7ra))2 
t  sin(7ra) 

s2  +  t2  +  2st  cos(7TQ;) 


The  proof  for  the  case  —  1  <  s  <  0  is  essentially  identical. 


(24) 


(25) 
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Corollary  2.6  Identities  (24)  and  (25)  remain  valid  after  any  rotation  or  translation 
of  the  curve  7  in  R2 . 

Corollary  2.7  When  the  curve  7  is  a  straight  line,  t/A  ^  ^(7 (t))  =  0  for  all  —  1  <  s  <  1 
and  —  1  <  t  <  1. 

2.3  Several  Classical  Analytical  Facts 

In  this  section  we  list  several  classical  analytical  facts.  They  can  be  found  in,  for  exam¬ 
ple,  [17]  and  [2], 

The  following  theorem  describes  a  property  of  the  zeros  of  analytic  functions. 

Theorem  2.8  If  f  is  a  nonzero  analytic  function  on  a  domain  Q  C  C,  then  the  zeros 
of  f  have  no  accumulation  point  in  Q. 

Corollary  2.9  (Analytic  continuation)  Suppose  that  f  and  g  are  both  analytic  func¬ 
tions  on  a  domain  O  C  t.  Suppose  further  that  f  and  g  are  equal  on  a  set  a  points  in 
fl  that  has  an  accumulation  point  in  11.  Then  f  and  g  are  equal  on  all  of  IT. 

The  following  classical  theorem  provides  a  test  for  the  convergence  of  an  infinite 
series. 

Theorem  2.10  (Dirichlet’s  test)  Suppose  a±,  <22, ...  is  a  sequence  of  real  numbers  such 
that 


^  ®ra+ 1  (28) 

for  each  positive  integer  n,  and 

lim  an  =  0.  (29) 

71— »•  OO 

Suppose  further  that  b±,  b2, . . .  is  a  sequence  of  complex  numbers  such  that,  for  some  real 
constant  M, 


N 


n=  1 


<  M, 


for  each  positive  integer  N.  Then 


(30) 


OO 

y,  anbn  <  00.  (31) 

n= 1 
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The  following  theorem  relates  a  limit  of  a  power  series  to  the  sum  of  its  coefficients. 

Theorem  2.11  (Abel’s  theorem)  Suppose  that  a0,  cp,  a2,  ■  ■ .  is  a  sequence  of  real  num¬ 
bers  such  that 

OO 

an£n  <  oo,  (32) 

n= 0 

for  all  —  1  <  x  <  1.  Suppose  further  that 

OO 

an  <  oo.  (33) 

n= 0 


T/ien 


lim  > 

i->i 

£<1  n.=0 


anxn  = 


E' 

n=0 


(34) 


3  Analytical  Apparatus 

The  elementary  theorems  3.4  and  3.5  in  this  section  are  the  primary  analytical  tools  of 
this  paper. 

The  following  theorem  provides  the  value  of  a  certain  integral.  It  is  found  in,  for  ex¬ 
ample,  [8],  Section  3.252,  formula  12.  For  completeness,  a  proof  is  provided  in  Appendix 
A. 


Theorem  3.1  Suppose  that  —l<p<l  and  0  <  a  <  2  are  real  numbers.  Then 


aAsm(7ra;)  ,  sm(/t7r(l  —  a)) 

- 7 - - - -  ax  =  7T t — r — ^ 

2 ax  cos(7tq;)  +  x2  sin(/t7r) 


(35) 


for  all  a  >  0. 


The  following  lemma  gives  the  Taylor  series  of  a  certain  rational  function. 


Lemma  3.2  Suppose  that  —  1  <  p  <  1  and  x  are  real  numbers.  Then 


p  sin(x) 

1  —  2pcos(x)  +  p2 


OO 

sin  (nx). 

n= 1 


(36) 
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Proof.  Let  —  1  <  p  <  1  and  x  be  real  numbers.  Then 


X>” 

n=  1 

=  Im  ( 


sm  rur  = 


lm[y^pneinx  }  =  Im 


rc=0 


1  —  pel 


1  —  pe 


\1  —  2pcos(a;)  +  p2 
p  sin(x) 

1  —  2pcos(a;)  +  p2 


(37) 


The  following  lemma  evaluates  the  integral  in  (35)  when  it  is  taken  from  0  to  1  instead 
of  from  0  to  oo. 

Lemma  3.3  Suppose  that  —  1  <  p  <  1  and  0  <  a  <  2  are  real  numbers.  Then 

[i ,  *"*»(**)  ;A=?ra.-sinMi^))+f  «-((* +d»«) (38) 

J o  a2  —  2ax  cos(7ra:)  +  x2  sm(/i7r) 

for  all  0  <  a  <  1. 

Proof.  Suppose  that  0  <  a  <  1.  Clearly, 

x^  sin(7ro;) 


k= 0 


p  —  k  —  1 


/,  a?  —  2 ax  cos(7r o)  +  x2 


dx  = 


(■OO 


^jsunyraj 


li  a  (f)2 -2(f)cos(vro)  + 1 


da;. 


(39) 


Since  -  <  1  for  all  x  >  1,  by  Lemma  3.2  we  observe  that 


a;M  sin(7TQ() 


L  a2  —  2aa;  cos(7ro;)  +  a;2 


da;  = 


OO 


E 

n=l 


an  sin(n7ro;) 


a;'1 


dx. 


(40) 


Interchanging  the  order  of  integration  and  summation,  we  further  observe  that 


a;M  sin(7ra) 


1 1  a2  —  2ax  cos(na)  +  x2 
an_1  sm(mra) 


dx  = 


°°  /‘°°  an  1  sm(nna) 


n= 1 


X 


n—fL+l 


dx 


p  —  n 


-E 

n= 1 

sin((A;  +  l)7ra;)^fe 


k= o 


p  —  k  —  1 


(41) 
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Combining  (35)  and  (41),  we  find  that 


x^  sin(7rad 


a2  —  2 ax  cos(7ra)  +  x 2 


,  sinfmRl  —  a)) 

dx  =  7T  a^1 K—, — +  Yl 

k= 0 


sin  (Jk  +  l)7ra) 


sin(/i7r) 


T 


(42) 


for  all  0  <  a  <  1. 


The  following  two  theorems  are  the  primary  analytical  tools  of  this  paper. 

A  simple  analytic  continuation  argument  shows  that  identity  (38)  in  lemma  3.3  is 
also  true  for  all  complex  /i  such  that  Re(/i)  >  —1  and  //  7^  1,  2,  3, . . ..  This  observation 
is  summarized  by  the  following  theorem. 


Theorem  3.4  Suppose  that  0  <  a  <  2  is  a  real  number  and  p  zs  complex,  so  that 
Re  fa  >  1  and  /z  7^  1,2,3, . . ..  Then 


oT  sin(7ro;) 


/0  a2  —  2ax  cos(7ra)  +  x2 


dx  =  n 


_1sin(/i7r(l  —  a))  sin((/c  +  l)7ra;) 


sin(/T7r) 


+  £ 

k= 0 


p  —  k  —  1 


(43) 


for  all  0  <  a  <  1. 


Proof.  Suppose  that  0  <  a  <  1.  We  observe  that  the  right  and  left  hand  sides  of 
identity  (38)  are  both  analytic  functions  of  //,  for  all  p  such  that  Re(/z)  >  —1  and 
/i  /  1,2,3,. . ..  Therefore,  by  analytic  continuation  (Theorem  2.9),  it  follows  that  iden¬ 
tity  (38)  holds  for  all  complex  p  such  that  Re(p)  >  —  1  and  p  7^  1,  2,  3, . . ..  ■ 


The  following  theorem  extends  Theorem  3.4  to  the  case  when  p  is  a  positive  integer. 
We  prove  it  by  repeated  application  of  L’Hopitahs  rule. 


Theorem  3.5  Suppose  that  0  <  a  <  2  is  a  real  number  mid  that  m  =  1,  2,  3, 
ri 


Then 


a:  sin  to 


dx 


Jo  °2  ~  Zax  cos(7ra)  +  x2 

Slid  \  \  rv  I  )7 TQf  ) 

=  am_17r(l  —  a)  cosOmna)  —  am_1  log(a)  sin(m7ro;)  +  — - - - — —ak, 

^ '  rn  —  k  —  1 


k>  0 

ky^m—1 


(44) 


for  all  0  <  a  <  1. 
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Proof.  Suppose  that  m  is  a  positive  integer  and  0  <  a  <  1.  We  observe  that  the  left 
hand  side  of  identity  (43)  is  analytic  in  /i,  for  all  p  such  that  Re(/i)  >  —1,  including  the 
points  p  =  1,  2, 3, . . ..  Therefore, 


dx  =  lim 


x M  sin(7ra;) 


o  a2  —  2 ax  cos(7ra;)  +  x2  J0  a2  —  2aa;  cos(7ro;)  +  x2 


dx.  (45) 


By  Theorem  3.4, 

rl 


'o 


ix  =  lim  TO.-sm(^(i-°))  +  £  »"((niH,» 


a2  —  2ax  cos(7tq;)  +  x2  m ->m  sin(/i7r)  '  fi  —  k  —  1 

sin((fc  +  l)vra)  |  _ ^sin^l  -  a))  |  sin(mTra)  nm_x 

/  7  Cl  I  ill'll  7 T CL  m  /  ^  | 


fc=0 

ky^m—1 


m  —  k  —  1 


fi—>m 


sm(pn) 


p  —  m 


V-  sin((fc  +  l)vra)  fc  p(/i) 

>  - 7 - t — a  +  Inn  — — , 

m  —  k  —  1  q{n ) 


k= 0 
k^m—1 


where  p  and  g  are  defined  by  the  formulas 

p(p)  =  7T a/i_1(/i  —  m)  sin(/x7r(l  —  a))  +  a™^1  sin(/i7r)  sin(m7ra:), 

and 

g(/i)  =  (p  —  m)  sin(/i7r). 

Clearly, 


lim  p(p)  =  lim  q(p)  =  0, 


so  we  will  use  L’Hopital’s  rule  to  determine  lim^™  We  observe  that 

p\p)  =  7rlog (a)aM_1(/x  —  m)  sin(/i7r(l  —  a))  +  7r aM_1  sin(/i7r(l  —  a)) 

+  7t2(1  —  a)aM_1(/i  —  m)  cos(yU7r(l  —  a))  +  7ram_1  cos(yU7r)  sin(m7ra),  (50) 


and 


Since 


q\p)  =  sin(/i7r)  +  (p  —  m)7i  cos(pn) . 


(46) 


(47) 


(48) 


(49) 


(51) 


lim  p\p)  =  lim  q\p)  =  0, 

H^m  fi—>m 


(52) 
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we  will  use  L’Hopital’s  rule  again.  We  observe  that 

p"{n)  =  7r(log(a))  W‘_1(yU  —  m)  sin(/i7r(l  —  a)  +  27rlog (a)aM_1  sin(/i7r(l  —  a)) 

+  27t2(1  —  a)  log (a)aM-1(/i  —  m)  cos(/i7r(l  —  a))  +  7t2(1  —  a)*/1-1  cos(yU7r(l  —  a)) 

—  7t3(1  —  a)2ah^1(/u  —  m)  sin(yU7r(l  —  a))  —  n2am~l  sin(/i7r)  sin(m7ra),  (53) 

and 

q'\p)  =  27rcos(/i7r)  —  (fi  —  m)^2  sin(/r7r).  (54) 


Since 


lim  p"(p)  —  2(— l)mam  17r2(l  —  a)  cos(rmra)  —  2(— l)mam  17r  log(a)  sin(m7ra),  (55) 

/i—^m 

and 


lim  g"(/j)  =  2vr(-l)m, 

H^m 


(56) 


it  follows  that 
/■i 


—  dx  —  ^  A 
/0  a 2  —  2ax  cos(7TQ;)  +  a:2  m  —  k  —  1 

k^m— 1 

v-  sin((/c  + l)7Ttt)  fc  p"(p) 

=  2^  - v - - — a  +  lim 


sin((fc  +  l)7ra)afc  +  ^  p{p) 

q^fl) 


k>  0 
k^m—l 


m  —  k  —  1 


M -m  q"(p) 

sm  ( [  K  |  j_  )7vcx ) 

=  am_17r(l  —  a)  cosfmna )  —  am_1  log(o)  sin(m7ro;)  +  - - - - — < 

^ '  rn  —  /d  —  1 


k>  0 
k^m—1 


(57) 


The  following  lemma  states  that  a  certain  series  converges. 

Lemma  3.6  Suppose  that  m  is  a  positive  integer  and  0  <  a  <  2  is  a  real  number. 

OO  .  /  \ 


£ 

n=  1 


sm  [nna , 
m  —  na 


<  oo. 


Then 

(58) 
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Proof.  We  observe  that 


> 


>0, 


net  —  m  (■ n  +  l)a  —  m 

for  all  positive  integers  n  such  that  n  >  m/a.  Moreover, 


lim  - =  0. 

n—too  na  —  m 


We  also  observe  that,  for  any  positive  integer  N, 

\  —  e  i(N+l)a 


N 

E 

n= 1 


sin  [nna 


< 


N 

E 

71= 0 


< 


0i(N+l)a/2  -i(N+l)a/2  _  i(N+l)a/2 


oia/2 


0—ia/2  _  pia/2 


_  cos((N  +  l)a/2) 
cos(o;/2) 

2 

<  - . 

|  cos(a/2)| 

Hence,  (58)  follows  by  Dirichlet’s  test  (Theorem  2.10). 


(59) 


(60) 


(61) 


The  following  theorem  states  that  a  certain  Taylor  series  converges  and  is  bounded 
on  the  interval  [0, 1]. 


Theorem  3.7  Suppose  that  m  and  k  are  positive  integers  and  0  <  a  <  2  is  a  real 
number.  Suppose  further  that  (p  is  a  function  [0, 1]  — *  R  defined  by  the  formula 


sin  {nna)tn_k 
m  —  na 


(62) 


for  all  0  <  t  <  1.  Then  f>  is  well  defined  and  bounded  on  the  interval  [0, 1]. 


Proof.  We  observe  that 


sin(7rncK) 

sin(7r(m  —  na)) 

-  •7T 

sin(7r(m  —  na)) 

m  —  na 

m  —  na 

-  /I 

7r(m  —  na) 

for  all  positive  integers  n.  Therefore, 


E 


sin  (t ma)tn_k 
m  —  na 


<  oo, 


(63) 


(64) 
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for  all  0  <  t  <  1.  By  Lemma  3.6, 


E 

n=k 


sm(7mm 


m  —  na 


<  oo, 


so  0  is  well  defined  on  [0, 1].  Furthermore,  by  Abel’s  theorem  (Theorem  2.11), 

H m  sm(TOa)r-t=^ 

t — ^-1  t  ^  m  —  r)  rv  t  J 


t->  i  z — '  m  —  na 

t<  1  n=k 


n=k 


m  —  na 


(65) 


(66) 


so  0  is  continuous  on  the  interval  [0, 1].  Therefore,  0  is  bounded  on  [0, 1]. 


The  following  theorem  states  that  a  certain  matrix  is  nonsingular. 


Theorem  3.8  Suppose  that  0  <  a  <  2  is  a  real  number.  Suppose  further  that  A(a )  is 
an  n  x  n  matrix  defined  via  the  formula 


sm  7rm 


a  ■ 


2j  —  1  —  ai 

sin(7rm) 

(2 -a) 


2 j  -  (2  -  a)i 


if  i  is  odd, 
if  i  is  even, 


(67) 


where  1  <  i,  j  <  n  are  integers.  Then  A(a )  is  nonsingular  for  all  but  a  finite  number  of 
0  <  a  <  2. 


Proof.  We  observe  that  the  functions 

sin(7rm) 

a  ■ 


(2 


2j  —  1  —  ai  ’ 

.  sin(7rm) 
-  a) 


2 j  -  (2  -  a) 


a)  i 


(68) 

(69) 


are  entire  functions  of  a,  where  1  <  i,j  <  n  are  integers.  Therefore,  det(A(a))  is  an 
entire  function  of  a.  We  also  observe  that 


A(l)  =  t rl,  (70) 

where  I  is  the  identity  matrix,  from  which  it  follows  that 

det(A(l))  =  7T.  (71) 

Since  the  interval  [0,2]  is  compact,  it  follows  from  Theorem  2.8  that  det(A(a))  is  equal 
to  0  at  no  more  than  a  finite  number  of  points  in  [0,2],  Hence,  A(a)  is  nonsingular  for 
all  but  a  finite  number  of  0  <  a  <  2. 
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4  Analysis  of  the  Integral  Equation:  the  Neumann 
Case 

Suppose  that  the  curve  7:  [—1, 1]  — >  R2  is  a  wedge  defined  by  (23)  with  interior  angle 
na,  where  0  <  a  <  2  (see  Figure  3).  Let  g  be  a  function  in  L2([—  1, 1]),  and  suppose  that 
p  G  L2([— 1, 1])  solves  the  equation 

-np(s)  +  J  ^{s)Hs){l(t))p(t)dt  =  g{s),  (72) 

for  all  s  G  [—1,1]. 

In  this  section,  we  will  analyze  this  boundary  integral  equation,  which  is  well-posed 
even  though  the  curve  7  is  open  (see  Observation  2.2).  In  Section  4.1  we  investigate  the 
behavior  of  (72)  for  functions  p  G  L2([— 1, 1])  of  the  forms 

Pit)  =  i*r\ 

pit)  =  sgn(t)IC~\ 

where  p  >  |  is  a  real  number  and 

(  —  1  if  x  <  0, 

sgn(x)  =  <  0  if  x  =  0,  (75) 

i  1  if  x  >  0, 

for  all  real  x.  If  identities  (24)  and  (25)  are  substituted  into  (72)  and  p  has  the  forms  (73) 
and  (74),  then  for  most  values  of  p  the  resulting  g  is  singular.  In  Section  4.2,  we  observe 
that  for  certain  p,  the  function  g  is  smooth.  In  Section  4.3,  we  fix  g  and  view  (72)  as 
an  integral  equation  in  p.  We  then  observe  that  for  certain  classes  of  functions  g,  the 
solution  p  is  representable  by  a  series  of  functions  of  the  forms  (73)  and  (74). 

4.1  Integral  Equations  Near  a  Corner 

The  following  lemma  uses  a  symmetry  argument  to  reduce  (72)  from  an  integral  equation 
on  the  interval  [—1, 1]  to  two  independent  integral  equations  on  the  interval  [0, 1]. 

Theorem  4.1  Suppose  that  p  is  a  function  in  L2([— 1,1])  and  that  g  G  L2([— 1,1])  is 
given  by  (72).  Suppose  further  that  even  functions  ge,  pe  G  L2([— 1, 1])  are  defined  via  the 
formulas 


1 

9eis)  =  2  (g(s) 

+  9(~s))i 

(76) 

1 

Peis)  =  ^  (d(S) 

+  p(-s)). 

(77) 

(73) 

(74) 
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Then 


9e(s) 


-TTpe(s) 


isin(7ra:)  .  .  , 

-^—72 — 71 - 7 — SPe^  dt 

s-  +  t2  —  1st  cos(7TQ;J 


(78) 


for  all  0  <  s  <  1. 

Likewise,  suppose  that  odd  functions  gQ,  pQ  €  L2([— 1, 1])  are  defined  via  the  formulas 


1 

9o(s)  =  ~(g(s)  -g(-s)), 

(79) 

Po(s)  =  \{p{s)  -  p(-s)). 

(80) 

Then 


9o(s) 


-npo(s) 


tsin(7ra)  .  .  , 

2  ~  - 7-, - 7 - 7P0 (t)  dt, 

s2  +  t2  —  2 st  COS(7TQ;) 


(81) 


for  all  0  <  s  <  1. 


Proof.  By  Lemma  2.3, 

.  .  .  f1  tsm(na)  .  .  , 

9{s)  =  ‘^(S)  -  f  s2  +  t2  +  2st  cos(Tra)  9  ^ 

for  all  —  1  <  s  <  0,  and 


0(s) 


-7rp(s)  + 


t  sin(7TQ() 


f_i  s2  +  t2  +  2st  cos(7rcx) 


p(t)  dt, 


for  all  0  <  s  <  1.  Therefore, 

.  .  .  f1  tsm(na )  .  . 

g(-s)  =  -M-*)  -  f  s2  +  (2_2s(cos(TO)l>W*, 

for  all  0  <  s  <  1,  and 

.  .  .  f1  tsin(na)  ,  . 

9(S)  =  -  l  s2  +  t2-2StcOS(„a)f‘{-t)  dt’ 

for  all  0  <  s  <  1. 

Adding  equations  (84)  and  (85),  we  observe  that 


9e(s) 


-npe(s) 


tsin(7ro;)  .  .  , 

~V~7 - 71 - 7 - ZPe(t)  dt 

S2  +t2  —  1st  COS(7rCK) 


(82) 


(83) 


(84) 


(85) 


(86) 
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(87) 


for  all  0  <  s  <  1. 

Likewise,  subtracting  equation  (84)  from  equation  (85),  we  observe  that 


9o(s)  =  -7r  pa(s) 


t  sin(7ra) 


/0  s2  +  t2  —  2 st  cos(7tq;) 


p0(t)  dt, 


for  all  0  <  s  <  1. 


4.2  The  Singularities  in  the  Solution  of  Equation  (72) 

In  this  section  we  observe  that  for  certain  functions  p,  the  functions  ge  and  g0  defined 
by  (76)  and  (79)  are  representable  by  convergent  Taylor  series  on  [0, 1]. 


4.2.1  The  Even  Case 

Suppose  that  p  G  L2([— 1,1])  is  an  even  function,  and  suppose  that  g  G  L2([— 1,1])  is 
defined  by  (72).  By  Theorem  4.1,  g  is  also  even  and 


9(s) 


-TTp(s) 


t  sinf-m) 


s2  +  t2  —  2 st  cos(7TQ;) 


p{t)  dt, 


(88) 


for  all  0  <  s  <  1. 

Suppose  further  that  p(t)  =  t for  all  0  <  t  <  1.  The  following  theorem  shows  that 
for  certain  values  of  //,  the  function  g  in  (88)  is  representable  by  a  convergent  Taylor 
series  on  the  interval  [0, 1]. 


Theorem  4.2  Suppose  that  0  <  a  <  2  is  a  real  number  and  n  is  a  positive  integer. 
Then 


zn—  l _ -i 

ITS  “  + 


zn  _ -I 

7 rs2-“  + 


t  sin(7ro;) 


!0  s2  +  f2  —  2st  cos(7tq;) 
t  sin(7ra;) 


I o  s2  T  t2  —  2 st  cos(7tq;) 


t2- 


-l 


t-i 


dt 


dt  =  (2 


m=  1 


sin(mTra)  ^m_1 
2n  —  1  —  am 


V  -  sin(-mTm)  x 

2n  —  (2  —  a)m 

m= 1  x  7 


(89) 

(90) 


for  all  0  <  s  <  1. 
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Proof.  Suppose  that  is  not  an  integer.  Substituting  p  =  into  (43),  we  observe 
that 


zrt—  1 _ -| 

ns  “  + 


2n  — 1  .  ,  \ 

s-  sm(7ra;) 


dt 


=  7TS 


=  7TS 


=  7TS 


t-1 


t-1 


^-1 


o  s2  —  2st  cos(7TQ;)  +  t2 

2n=i_iSin(^  •  7r(l  -  a)) 


+  ns 


+  ns 


sm(^-n) 


sinf^i  -n) 


sin ((fc  +  1)7TQ;) 

+  Z^  2n=l_fc_1 

Esin((h,  - 


fc=0  « 


2n-i_1sin(^^  •  7T  —  (2n  —  l)7r)  sin((/c  +  l)7ro;)  k 


k= o 


2n  —  1  —  a(k  +  1) 


—  ns 


OC 


sin(m7ra) 


0ra— 1 


^  2n  —  1  —  am 

m=  1 


=  a 


sin  (run a) 

2^  9  n  _  i  _ 


1 


m=l 


2n  —  1  —  am 


(91) 


for  all  s  >  0. 

Now  suppose  that  -  is  an  integer.  We  observe  that  there  is  a  neighborhood  V  of 
a  such  that  is  not  an  integer  on  V\  {«}.  Clearly,  (91)  is  a  bounded  and  analytic  on 
V\ {ck}.  Therefore,  identity  (91)  is  extended  to  this  case  by  an  application  of  L’Hopital’s 
rule. 


4.2.2  The  Odd  Case 

Suppose  that  p  G  L2([— 1,1])  is  an  odd  function,  and  suppose  that  g  G  L2([— 1,1])  is 
defined  by  (72).  By  Theorem  4.1,  g  is  also  odd  and 


9{s) 


-n  p(s) 


t  sin(7ro;) 


s2  +  t2  —  2stcos(na) 


pit )  dt, 


(92) 


for  all  0  <  s  <  1. 

Suppose  further  that  p(t)  =  t M_1  for  all  0  <  t  <  1.  The  following  theorem  shows  that 
for  certain  values  of  //,  the  function  g  in  (92)  is  representable  by  a  convergent  Taylor 
series  on  the  interval  [0, 1]. 


Theorem  4.3  Suppose  that  0  <  a  <  2  is  a  real  number  and  n  is  a  positive  integer. 
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Then 


—  7 rs 


l--i 


t  sin(7ro;) 


2n—l  1 

—  7 rs  2-“  + 


/ o  s2  +  t2  —  2st  cos(7ra) 
t  sin(7ra:) 


2 n  !  ,  sin(m7ra)  „  . 

t^~1dt  =  a>  — ^ - -sm~\ 


^  2 n  —  am 

m=  1 


(93) 


*0  s2  +  t2  —  2 si  cos(7tq;) 


t 


t-i 


dt=(2-a)Y - 

v  ^  9rj  -  1  -  9.  - 


Tm— 1 


771=1 


2n 


(2 


arm 


(94) 


for  all  0  <  s  <  1. 


4.3  Series  Representation  of  the  Solution  of  Equation  (72) 

Suppose  that  g  is  a  function  in  L2([— 1, 1])  such  that  ge  and  ga  defined  by  (76)  and  (79)  are 
representable  by  convergent  Taylor  series  on  [0, 1].  Suppose  further  that  p  G  L2([— 1, 1]) 
satisfies  equation  (72).  In  this  section  we  observe  that  pe  and  pa  defined  by  (77)  and  (80) 
are  representable  by  certain  series  of  singular  powers  of  x  G  [0, 1]. 


4.3.1  The  Even  Case 


Suppose  that  g  G  L2([— 1,1])  is  an  even  function,  and  suppose  that  p  G  L2([— 1,1]) 
satisfies  equation  (72).  By  Theorem  4.1,  p  is  also  even  and 


-  Ms) 


t  sin(7ra) 

s2  +  t2  —  2 st  cos(7tq;) 


p(t)dt  =g(s), 


(95) 


for  all  0  <  s  <  1,  where  0  <  a  <  2. 

Let  |" x]  denote  the  smallest  integer  n  such  that  n  >  x,  and  let  |ycj  denote  the  largest 
integer  n  such  that  n  <  x,  for  all  real  x.  The  following  theorem  shows  that  if  the  g  is 
representable  by  a  convergent  Taylor  series  on  [0, 1],  then  for  any  positive  integer  n  there 
exist  unique  real  numbers  b\,  62,  •  •  • ,  bn  such  that  the  function 


\n/2\  |n/2j 

pit)  =  y  b2i-it2jL^~i  +  y 

7=1  7=1 


(96) 


where  0  <  t  <  1,  solves  equation  (95)  to  within  an  error  0(tn). 


Theorem  4.4  Suppose  that  n  is  a  positive  integer  and  ci,c2, . . . ,  cn  are  real  numbers. 
Suppose  further  that  g:  [0, 1]  — >■  R  is  defined  by  the  formula 

n 

g(t)  =  Yc^  ^ 

7=1 
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for  all  0  <  t  <  1.  Then,  for  all  but  a  finite  number  of  0  <  a  <  2,  there  exist  unique  real 
numbers  bi,  b2,  . . . ,  bn  such  that 


\n/2\  |n/2j 

p{t)  =  22  +  22 h<2itl 

i= 1  2=1 


-1 


(98) 


t  sin(7ra:) 


for  all  0  <  t  <  1,  and 

^  ^  J0  s2  +  t2  —  2st  cos(7ra)y 

for  all  0  <  s  <  l,  where  (f:  [0, 1]  — >  R  is  a  bounded  function  representable  by  a  convergent 
Taylor  series  of  the  form 


-pit)  dt  =  g(s)  +  sn<j>(s), 


(99) 


(fit)  =  22 dit 1  \  (100) 

2=1 

for  all  0  <  t  <  1,  where  di,  d2,...  are  real  numbers. 

Proof.  By  Theorem  3.8,  the  n  x  n  matrix  A{a)  defined  by  (67)  is  nonsingular  for  all 
but  a  finite  number  of  0  <  a  <  2.  Whenever  A{a)  is  nonsingular,  there  exist  unique  real 
numbers  b\,  b2, . . . ,  bn  such  that 


^  ^  Ai®)ijbj  Q, 

j- 1 


(101) 


for  every  i  —  1,  2, . . . ,  n.  Suppose  that  p :  [0, 1]  — *  R  is  defined  by  (98).  By  Theorem  4.2, 
f1  t  sin(7ro;) 


vrp(s) 


/0  s2  +  t2  —  2 st  COs(7TCk) 


pit)  dt 


fn/21 


sin  ra 


Ln/2j 


= E  ^  E  2;_  i_m; 

j= 1  *=1 

fn/21  n  .  /  *  \ 

E,  sm(7rcm 

^  E 

j=i  i=i  J 


1  -  (2  -  a)  22  21  W 


sm  7rm 


_2 —  1 


-  J  “j'  2j  —  (2  —  a)i 


3 

Ln/2j 


1  -  (2  -  a)  ^  b2jJ2  w~. 


sm  7rm 


-  ^^2j-(2-a)< 


V-i+s^(s), 


(102) 


for  all  0  <  s  <  1,  where  <f :  [0, 1]  — »  R  is  defined  by  the  formula 

fn/21  co  .  /  .\  Ln/2J  oo 


0(t)  =  -a  62j-i 


sm  7rcm 


j=l  2=n+l 


2  j  —  1  —  m 


-  (2  -  a)  E  ^  E 


sm  nrm 


#2— 1— n 


j= 1  2=n+l 


2j  -  (2  -  a)/ 


(103) 
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for  all  0  <  t  <  1.  By  Theorem  3.7,  (j)  is  bounded  on  [0, 1].  By  interchanging  the  order  of 
summation  in  (102),  we  observe  that 


vrp(s) 


t  sin(7ro;) 


J o  s2  +  t2  —  2 st  cos(7TQ;) 

fn/21 


pit )  dtt 


E  HE 


sm^ai 
— (  2j'  -  1  -  od 


Ln/2j 


T&2j-i  ~  (2  -  a)  T 


sm  [Tim) 


-b2j)  sl  1  +  sncj)(s ) 

^  X  ' 


2 j  -  (2  -  a)t 

+  s”‘#,W  =  ECiS*_1  +  S’V(S)  =  sW  +  s’VMi  (104) 


2=1  J  =  1 

for  all  0  <  s  <  1. 


2=1 


4.3.2  The  Odd  Case 

Suppose  that  g  G  L2([—  1, 1])  is  an  odd  function,  and  suppose  that  p  G  L2([— 1, 1])  satisfies 
equation  (72).  By  Theorem  4.1,  p  is  also  odd  and 


-  7T p(s) 


t  sin(7ra) 

s2  +  t2  —  2 st  cos(7tq;) 


p(t)dt  =  g(s), 


(105) 


for  all  0  <  s  <  1,  where  0  <  a  <  2. 

Let  denote  the  smallest  integer  n  such  that  n  >  x,  and  let  denote  the  largest 
integer  n  such  that  n  <  x,  for  all  real  x.  The  following  theorem  shows  that  if  the  g  is 
representable  by  a  convergent  Taylor  series  on  [0, 1],  then  for  any  positive  integer  n  there 
exist  unique  real  numbers  b\,  b2,  ■  ■  ■ ,  bn  such  that  the  function 


\n/ 21 

P(t) 


-d2-" 


2j— 1 


Ln/2J 

-1  +  E b 


2  U 


Zl  -1 


2=1 


2=1 


(106) 


where  0  <  t  <  1,  solves  equation  (105)  to  within  an  error  0(tn). 


Theorem  4.5  Suppose  that  n  is  a  positive  integer  and  ci,c2, . . . ,  cn  are  real  numbers. 
Suppose  further  that  g :  [0, 1]  — »  R  is  defined  by  the  formula 

n 

9(t)  =  ^2<ktt~1,  (107) 

2=1 
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for  all  0  <  t  <  1.  Then,  for  all  but  a  finite  number  of  0  <  a  <  2,  there  exist  unique  real 
numbers  bi,  b2, . . . ,  bn  so  that 

\n/2\  \n/2\ 

pit)  =  y  +  Y  62i<  “-1,  (108) 

i= 1  2—1 

for  all  0  <  t  <  l,  and 

+  l  S2  +  (2!l7cL(TO)',(t)  ^  =  J(S)  +  S>(S)'  (1M) 

for  all  0  <  s  <  1,  where  0:  [0, 1]  — »  R  is  a  bounded  function  representable  by  a  convergent 
Taylor  series  of  the  form 

OO 

0(*) = y:  djti~ir  (no) 

2=1 

for  all  0  <  t  <  1,  where  d\,  d2,  •  •  •  are  real  numbers. 


4.4  Summary  of  Results 

We  summarize  the  results  of  the  preceding  subsections  4.1,  4.2,  4.3  as  follows. 

Suppose  that  the  curve  7 :  [—1, 1]  — >  R2  is  a  wedge  defined  by  (23)  with  interior  angle 
na,  where  0  <  a  <  2  (see  Figure  3).  Let  g  G  L2([— 1,1]),  and  consider  the  boundary 
integral  equation 

-vrp(s)  +  J  ip^8)tV{s)('y(t))p(t)dt  =  g(s),  (111) 

for  all  s  G  [—1, 1],  where  p  G  L2([— 1, 1]). 

Suppose  that  the  even  and  odd  parts  of  g  are  each  representable  by  convergent  Taylor 
series  on  the  interval  [0, 1].  Then,  for  each  positive  integer  n,  there  exist  real  numbers 
61,  62, . . . ,  bn  and  <7,  c2, . . . ,  cn  such  that 


\n/2\  [n/2j 

p(t)  =  Y  +  Y  &2iSgn(t)|i| 

2=1  2=1 


K21 

+  Y  c2i-iSgn(t)|t|^_1 


Ln/2j 

+  Y 

2=1 


(112) 


for  all  —1  <  t  <  1,  solves  equation  (111)  to  within  an  error  0(tn).  Moreover,  the  even 
and  odds  parts  of  this  error  are  also  representable  by  convergent  Taylor  series  on  the 
interval  [0, 1]  (see  theorems  4.4  and  4.5). 
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Observation  4.1  Numerical  experiments  ( see  Section  6)  suggest  that,  for  a  certain  sub¬ 
class  of  functions  g,  stronger  versions  of  theorems  4-4  and  4-5  are  true.  Suppose  that  G 
is  a  harmonic  function  on  a  neighborhood  of  the  set  {x  G  R2  :  ||x||  <  1},  and  let 

f)G 

m  =  mm),  (H3) 

for  all  —  1  <  t  <  l,  where  v{t)  is  the  inward-pointing  unit  normal  vector  at  7 (t).  We 
conjecture  that  there  exist  infinite  sequences  of  real  numbers  bi,  62,  •  •  •  and  <7,  C2, . . .  such 
that 


p(t)  =  b2i~i 

i= 1 


2Iq1  1  +  J]&2iSgn(t)|t|«  1 

i= 1 


00 

+  J^C2i-iSgn(t)|t|^ 

2—1 


(114) 


2—  1 


is  well  defined  for  all  —  1  <  t  <  1,  and  (114)  solves  equation  (111). 


Observation  4.2  Numerical  experiments  (see  Section  6)  indicate  that  the  solution  to 
equation  (111)  is  representable  by  a  series  of  the  form  (112),  to  high  precision,  for  a  more 
general  class  of  curves  7.  More  specifically,  suppose  that  7:  [—1, 1]  — *  R2  is  a  wedge  in 
R2  with  smooth,  curved  sides,  with  a  corner  at  0  and  interior  angle  71a.  Suppose  further 
that  all  derivatives  of  7,  2nd  order  and  higher,  approach  zero  near  the  corner.  Then  the 
solution  is  representable  by  a  series  of  the  form  (112),  to  high  precision. 


5  Analysis  of  the  Integral  Equation:  the  Dirichlet 
Case 

Suppose  that  the  curve  7:  [—1, 1]  — >  R2  is  a  wedge  defined  by  (23)  with  interior  angle 
na,  where  0  <  a  <  2  (see  Figure  3).  Let  g  be  a  function  in  L2([— 1, 1]),  and  suppose  that 
p  G  L2([— 1, 1])  solves  the  equation 

9(s)  =  -Kp(s)  +  J  (115) 

for  all  .s  G  [—1,1]. 

In  this  section,  we  will  analyze  this  boundary  integral  equation,  which  is  well-posed 
even  though  the  curve  7  is  open  (see  Observation  2.2).  In  Section  5.1  we  investigate  the 
behavior  of  (115)  for  functions  p  G  L2([— 1, 1])  of  the  forms 

p(t)  =  111“. 

p(t)  =  Sgll(f)|«|'\ 


(116) 

(117) 


where  p,  >  1  is  a  real  number  and 

(  —  1  if  x  <  0, 

sgn(x)  =  <  0  if  x  =  0,  (118) 

[  1  if  x  >  0, 

for  all  real  x.  If  identities  (24)  and  (25)  are  substituted  into  (115)  and  p  has  the 
forms  (116)  and  (117),  then  for  most  values  of  //  the  resulting  g  is  singular.  In  Sec¬ 
tion  5.2,  we  observe  that  for  certain  p,  the  function  g  is  smooth.  In  Section  5.3,  we  fix 
g  and  view  (115)  as  an  integral  equation  in  p.  We  then  observe  that  for  certain  classes 
of  functions  g ,  the  solution  p  is  representable  by  a  series  of  functions  of  the  forms  (116) 
and  (117). 

The  proofs  of  the  theorems  in  this  section  are  essentially  identical  to  the  proofs  of 
the  corresponding  theorems  in  Section  4,  and  are  omitted. 


5.1  Integral  Equations  Near  a  Corner 

The  following  lemma  uses  a  symmetry  argument  to  reduce  (115)  from  an  integral  equation 
on  the  interval  [—1, 1]  to  two  independent  integral  equations  on  the  interval  [0, 1]. 

Theorem  5.1  Suppose  that  p  is  a  function  in  L2([— 1,1])  and  that  g  G  L2([— 1,1])  is 
given  by  (115).  Suppose  further  that  even  functions  ge,  pe  G  L2([— 1,1])  are  defined  via 
the  formulas 


9e{s) 

=  \(sM  +  «(-»)). 

(119) 

Pe(s) 

=  1(pW  +  p(-«)). 

(120) 

9e{s) 

fl  ssinf7ro;)  .  , 

(121) 

=  nMs)  (  s^  +  r-2stL{^{t)dt' 

for  all  0  <  s  <  1. 

Likewise,  suppose  that  odd  functions  ga,  p0  G  T2([— 1, 1])  are  defined  via  the  formulas 
ffo(s)  =  ^(g(s)-g(-s)),  (122) 

Po(s)  =  ^(p(s)  -  p(-s)).  (123) 

Then 

.  .  .  .  f1  ssin(7ro;)  .  .  ,  . 

J.M  =  +  f  s2  +  ^2s(cos(to)p»M  *.  (124) 

for  all  0  <  s  <  1. 
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5.2  The  Singularities  in  the  Solution  of  Equation  (115) 

In  this  section  we  observe  that  for  certain  functions  p,  the  functions  ge  and  gQ  defined 
by  (119)  and  (122)  are  representable  by  convergent  Taylor  series  on  [0, 1]. 


5.2.1  The  Even  Case 

Suppose  that  p  G  L2([— 1,1])  is  an  even  function,  and  suppose  that  g  G  L2([— 1,1])  is 
defined  by  (115).  By  Theorem  5.1,  g  is  also  even  and 


9{s) 


-7T p(s)  ~ 


ssm  to 


s2  +  t2  —  2 st  cos(7ra) 


p{t)  dt, 


(125) 


for  all  0  <  s  <  1. 

Suppose  further  that  p(t)  =  t M_1  for  all  0  <  t  <  1.  The  following  theorem  shows  that 
for  certain  values  of  p,  the  function  g  in  (125)  is  representable  by  a  convergent  Taylor 
series  on  the  interval  [0, 1]. 


Theorem  5.2  Suppose  that  0  <  a  <  2  is  a  real  number  and  n  is  a  positive  integer. 
Then 


ssm  to 


n 


2n  —  l 

7TS  “  + 


l0  s2  +  t2  —  2 st  COs(7TCk) 
r  1 


dt  =  (2  —  o)ti  — 


m=  1 
oo 


s  sin  m 


2  n 

7TS2~a  + 


I o  s2  +  t2  —  2 st  COs(7TQ;) 
s  sin(7ro;) 


2n  — 1  V — 

t  “  dt  =  a  y 


m  =  1 


sinymiux)  m 
■S  ) 
m 

sin(m7ro;) 

2  n  —  1  —  am 


l0  s2  +  t2  —  2 st  cos(7tq;) 


t2~a  dt  =  (2  —  a) 


OO 


sin  mM 


^  2 n  -  (2 

m— 1  x 


am 


(126) 

(127) 

(128) 


for  all  0  <  s  <  1. 


Proof.  Taking  the  limit  p,  — >  0  in  (43)  and  applying  L’Hopital’s  rule  once,  we  observe 
that 


sin  7ra 


/0  a2  —  2 ax  cos(7tck)  +  x2 


dx  —  (1  —  a)na  — 


OO 

sin 


k= 0 


[fk  +  l)7ra) 

k  +  1 


(129) 


for  all  0  <  a  <  1.  Clearly,  identity  (126)  follows. 

The  proofs  of  identities  (127)  and  (128)  are  essentially  identical  to  the  corresponding 
proofs  in  Theorem  4.2.  ■ 
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5.2.2  The  Odd  Case 


Suppose  that  p  G  L2([— 1,1])  is  an  odd  function,  and  suppose  that  g  G  L2([— 1,1])  is 
defined  by  (72).  By  Theorem  4.1,  g  is  also  odd  and 


9{s) 


-irp(s) 


t  sin(7ro;) 


s2  +  t2  —  2stcos(na) 


p{t)  dt , 


(130) 


for  all  0  <  s  <  1. 

Suppose  further  that  p(t)  =  t M_1  for  all  0  <  t  <  1.  The  following  theorem  shows  that 
for  certain  values  of  //,  the  function  g  in  (92)  is  representable  by  a  convergent  Taylor 
series  on  the  interval  [0, 1]. 


Theorem  5.3  Suppose  that  0  <  a  <  2  is  a  real  number  and  n  is  a  positive  integer. 
Then 


-  vr  + 


ITS  a  T 


ssm  to 


2re— 1 

ITS  2~“  + 


l0  s2  +  t2  —  2st  cos(7ra) 

f1  ssin(7ro;) 

/0  s2  +  t2  —  2 st  cos(7TQ;) 

s  sin(7ro;) 


dt  =  —air  — 


m=  1 
oo 


t  a  dt  =  a 


m= 1 


smyrmra ) 
m 

sin(m7ra) 
2  n  —  am 


l0  s2  +  t2  —  2stcos(na) 


1 dt  =  (2  —  a) 


OO 


sin  [mira) 


(131) 

(132) 


^  In  -  1  -  (2 

m=  1 


am 


(133) 


for  all  0  <  s  <  1. 


5.3  Series  Representation  of  the  Solution  of  Equation  (115) 

Suppose  that  g  is  a  function  in  L2([— 1, 1])  such  that  ge  and  gQ  defined  by  (76)  and  (79)  are 
representable  by  convergent  Taylor  series  on  [0,1].  Suppose  further  that  p  G  L2([— 1, 1]) 
satisfies  equation  (72).  In  this  section  we  observe  that  pe  and  pQ  defined  by  (77)  and  (80) 
are  representable  by  certain  series  of  singular  powers  of  x  G  [0, 1]. 


5.3.1  The  Even  Case 

Suppose  that  g  G  L2([— 1,1])  is  an  even  function,  and  suppose  that  p  G  L2([— 1,1]) 
satisfies  equation  (72).  By  Theorem  4.1,  p  is  also  even  and 


-  7T p(s) 


s  sin(7ra) 

s2  +  t2  —  2 st  cos(7ra) 


p{t)dt  =  g(s), 


(134) 


31 


for  all  0  <  s  <  1,  where  0  <  a  <  2. 

Let  |V|  denote  the  smallest  integer  n  such  that  n  >  x,  and  let  |yrj  denote  the  largest 
integer  n  such  that  n  <  x,  for  all  real  x.  The  following  theorem  shows  that  if  the  g  is 
representable  by  a  convergent  Taylor  series  on  [0, 1],  then  for  any  positive  integer  n  there 
exist  unique  real  numbers  b0,bi, ... ,  bn  such  that  the  function 

Tn/21  |rc/2j 

pit)  =  22  hi-it^  +  22  b2it^,  (135) 

i— 1 

where  0  <  t  <  1,  solves  equation  (134)  to  within  an  error  0(tn+1). 

Theorem  5.4  Suppose  that  n  is  a  positive  integer  and  c0,ci, . . .  ,cn  are  real  numbers. 
Suppose  further  that  g:  [0, 1]  — *  R  is  defined  by  the  formula 

n 

g(t)  =  '22cit\  (136) 

i= 0 

for  all  0  <  t  <  1.  Then,  for  all  but  a  finite  number  of  0  <  a  <  2,  there  exist  unique  real 
numbers  b0,  bi, . . . ,  bn  so  that 

\n/2\  _  In/ 2] 

p{t)  =  ^2  b2i-\t^  +  22  b2it^,  (137) 

2=1  2=0 

for  all  0  <  t  <  1,  and 

Tipis)  f  0tSmJ7rQ^  , - Tp(t)dt  =  gis)  +  sn+V(s),  (138) 

J o  sz  +  tz  —  1st  cos(7TQ;) 

for  all  0  <  s  <  1,  where  0:  [0, 1]  — »  R  is  a  bounded  function  representable  by  a  convergent 
Taylor  series  of  the  form 

OO 

(fit)  =  22  ditl’  (139) 

i= 0 

for  all  0  <  t  <  1,  where  d0 ,  di, . . .  are  real  numbers. 


5.3.2  The  Odd  Case 


Suppose  that  g  G  L2([— 1, 1])  is  an  odd  function,  and  suppose  that  p  G  L2([— 1, 1])  satisfies 
equation  (72).  By  Theorem  4.1,  p  is  also  odd  and 


-  7T pis) 


s  sin(7ra) 

s2  +  t2  —  2 st  cos(7ra) 


p{t)dt  =  gis), 


(140) 
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for  all  0  <  s  <  1,  where  0  <  a  <  2. 

Let  |V|  denote  the  smallest  integer  n  such  that  n  >  x,  and  let  |ycj  denote  the  largest 
integer  n  such  that  n  <  x,  for  all  real  x.  The  following  theorem  shows  that  if  the  g  is 
representable  by  a  convergent  Taylor  series  on  [0, 1],  then  for  any  positive  integer  n  there 
exist  unique  real  numbers  b0,b\, ... ,  bn  such  that  the  function 

\n/2\  _  \n/2\ 

p(t )  =  Y  +  Y  ,  (141) 

2=  1  2=0 

where  0  <  t  <  1,  solves  equation  (140)  to  within  an  error  0(tn+1). 

Theorem  5.5  Suppose  that  n  is  a  positive  integer  and  co,ci, . . . ,  cn  are  real  numbers. 
Suppose  further  that  g:  [0, 1]  — *  R  is  defined  by  the  formula 

n 

g{t)  =  Y^t\  (142) 

i= 0 

for  all  0  <  t  <  1.  Then,  for  all  but  a  finite  number  of  0  <  a  <  2,  there  exist  unique  real 
numbers  b0,bi, ... ,  bn  so  that 

rn/21  K2J 

pit )  =  Y  hi-it^  +  Y  ,  (143) 

2=1  2=0 

for  all  0  <  t  <  1,  and 

-  Ms)  + 1  ^  +  e-n2(rl(TO/(t)  dt  =  9(s)  +  i”+V(s)-  (144) 

for  all  0  <  s  <  l,  where  fi:  [0, 1]  — »  R  is  a  bounded  function  representable  by  a  convergent 
Taylor  series  of  the  form 

OO 

(fit)  =  Y  (145) 

i= 0 

for  all  0  <  t  <  1,  where  d0,  di, . . .  are  real  numbers. 

5.4  Summary  of  Results 

We  summarize  the  results  of  the  preceding  subsections  5.1,  5.2,  5.3  as  follows. 
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Suppose  that  the  curve  7 :  [—1, 1]  — »  R2  is  a  wedge  defined  by  (23)  with  interior  angle 
7TQ;,  where  0  <  a  <  2  (see  Figure  3).  Let  g  G  L2([— 1,1]),  and  consider  the  boundary 
integral  equation 

*p(s)  +  J  ^S),u(s)('y(t))p(t)dt  =  g{s),  (146) 

for  all  s  G  [—1, 1],  where  p  G  L2([— 1, 1]). 

Suppose  that  the  even  and  odd  parts  of  g  are  each  representable  by  convergent  Taylor 
series  on  the  interval  [0, 1].  Then,  for  each  positive  integer  n,  there  exist  real  numbers 
60,  bn  and  c0,  <7, . . . ,  cn  such  that 


\n/2\  \n/2\  fn/2]  \n/2\ 

P(t)  =  ^2  b2i-i\t\^1  +  ^2  &2iSgn(t)|t|«  +  ^2  c2i-isgn(t)|t|^  +  ^  c2i\t\ 


2  i 

2  —  a 


2=1 


2=0 


2=1 


2=0 


(147) 


for  all  —  1  <  t  <  1,  solves  equation  (146)  to  within  an  error  0(tn+1).  Moreover,  the  even 
and  odds  parts  of  this  error  are  also  representable  by  convergent  Taylor  series  on  the 
interval  [0,1]  (see  theorems  5.4  and  5.5). 


Observation  5.1  Numerical  experiments  ( see  Section  6)  suggest  that,  for  a  certain  sub¬ 
class  of  functions  g,  stronger  versions  of  theorems  5.f  and  5.5  are  true.  Suppose  that  G 
is  a  harmonic  function  on  a  neighborhood  of  the  set  {a:  G  R2  :  ||a;||  <  1},  and  let 

gif )  =  G(j(t)),  (148) 


for  all  —  1  <  t  <  1.  We  conjecture  that  there  exist  infinite  sequences  of  real  numbers 
b0,  b\, . . .  and  cq,  <7, . . .  such  that 


OO  OO  OO 

P(t)  =  +^2b2isgn(t)\t\  «  +  ^c2i-isgn(t)|t|^ 

2=1  2=0  2=1 


+  '22c2i\t\ 2-“, 

i= 0 


(149) 


is  well  defined  for  all  —  1  <  t  <  1,  and  (lj9)  solves  equation  (lj6). 


Observation  5.2  Numerical  experiments  (see  Section  6)  indicate  that  the  solution  to 
equation  (lf6)  is  representable  by  a  series  of  the  form  (lf7),  to  high  precision,  for  a  more 
general  class  of  curves  7.  More  specifically,  suppose  that  7:  [—1, 1]  — >  R2  is  a  wedge  in 
R2  with  smooth,  curved  sides,  with  a  corner  at  0  and  interior  angle  1 ra.  Suppose  further 
that  all  derivatives  of  7,  2nd  order  and  higher,  approach  zero  near  the  corner.  Then  the 
solution  is  representable  by  a  series  of  the  form  (147),  to  high  precision. 
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6  The  Algorithm 

To  solve  the  integral  equations  of  potential  theory  on  polygonal  domains,  we  use  the 
algorithm  described  in  [3] ;  however,  instead  of  discretizing  the  corner  singularities  using 
nested  quadratures,  we  use  the  representations  (112),  (147)  to  construct  purpose-made 
discretizations  (see,  for  example,  [18],  [16],  [24]).  A  detailed  description  of  this  part  of 
the  procedure  will  be  published  at  a  later  date. 

We  illustrate  the  performance  of  the  algorithm  with  several  numerical  examples. 
The  exterior  Neumann  problem  and  interior  Dirichlct  problem  were  solved  on  each  of 
the  domains  in  figures  4-9,  where  the  boundary  data  was  generated  by  a  unit  charge 
inside  the  region  in  the  Neumann  case,  and  outside  the  region  in  the  Dirichlct  case.  The 
numerical  solution  was  tested  by  comparing  the  computed  potential  to  the  true  potential 
at  five  arbitrary  points.  Tables  1  and  2  present  the  results.  The  following  quantities  are 
reported  for  each  problem: 

n  is  the  total  number  of  nodes; 

t  is  the  wall  clock  time  required  to  solve  the  problem; 

E  is  the  largest  error  in  absolute  value  measured  in  the  computed  potential; 
k  is  the  condition  number  of  the  linear  system  that  was  solved. 

Observation  6.1  Clearly,  the  curves  Ti  and  T2  are  not  polygons.  However,  all  deriva¬ 
tives  of  the  curve,  2nd  order  and  higher,  approach  zero  near  the  corner.  We  observe 
that  in  this  case,  the  singularities  in  the  solutions  of  the  boundary  integral  equations  are 
identical  to  those  in  the  polygonal  case. 

Observation  6.2  We  observe  that  if  the  boundary  values  are  produced  by  a  charge  inside 
the  domain  in  the  Neumann  case,  or  outside  the  domain  in  the  Dirichlet  case,  certain 
terms  in  the  representations  of  the  solutions  near  a  corner  vanish.  More  specifically,  in 
the  Neumann  case,  the  terms  ci,c2, . . .  in  (Ilf)  vanish,  while  in  the  Dirichlet  case,  the 
terms  b0,  bi, . . .  in  (If  9)  vanish. 

Observation  6.3  The  purpose  of  the  numerical  examples  in  this  section  is  to  illustrate 
the  accuracy  and  convergence  of  the  algorithms,  as  opposed  to  the  resulting  CPU  time. 
For  historical  reasons,  the  resulting  linear  systems  were  solved  using  the  minimal  residual 
method  (see,  for  example,  [7]).  Clearly,  the  CPU  times  in  tables  1  and  2  would  be 
significantly  improved  by  using  an  algorithm  appropriate  to  the  size  of  the  problems  being 
solved. 
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7  Extensions  and  Generalizations 


7.1  The  Helmholtz  Equation  on  Polygonal  Domains 

In  this  report,  we  investigate  the  solution  of  boundary  value  problems  for  Laplace’s  equa¬ 
tion  on  polygonal  domains;  we  observe  that,  when  the  problems  are  formulated  as  the 
boundary  integral  equations  of  classical  potential  theory,  the  solutions  are  representable 
by  series  of  elementary  functions.  A  similar  analysis  applies  to  the  Helmholtz  equation 
on  polygonal  domains.  More  specifically,  if  boundary  value  problems  for  the  Helmholtz 
equation  on  polygonal  domains  are  formulated  as  the  boundary  integral  equations  of 
classical  potential  theory,  the  solutions  are  representable  by  series  of  appropriately  se¬ 
lected  Bessel  functions.  A  paper  containing  this  analysis  (together  with  the  requisite 
numerical  apparatus)  is  in  preparation. 

7.2  Curved  Boundaries  with  Corners 

While  this  report  only  deals  with  the  solution  of  Laplace’s  equation  on  domains  with 
polygonal  boundaries,  a  similar  analysis  applies  to  the  case  of  curved  boundaries  with 
corners.  More  specifically,  if  the  boundary  is  smooth  except  at  corners,  the  solutions 
to  the  associated  boundary  integral  equations  of  classical  potential  theory  are  also  rep¬ 
resentable  by  series  of  elementary  functions.  This  analysis,  along  with  the  requisite 
numerical  apparatus,  will  be  described  in  a  forthcoming  paper. 

7.3  Generalization  to  Three  Dimensions 

The  generalization  of  the  apparatus  of  this  report  to  three  dimensions  is  fairly  straight¬ 
forward,  but  the  detailed  analysis  has  not  been  carried  out.  This  line  of  research  is  being 
vigorously  pursued. 

7.4  Robin  and  Mixed  Boundary  Conditions 

This  report  deals  with  the  solution  of  Laplace’s  equation  on  polygonal  domains  with 
either  Dirichlet  or  Neumann  boundary  conditions.  There  are  two  additional  boundary 
conditions  that  have  not  yet  been  analyzed  in  detail:  the  Robin  condition,  which  specifies 
a  linear  combination  of  the  values  of  the  solution  and  the  values  of  its  derivative  on 
the  boundary;  and  the  mixed  boundary  condition,  which  specifies  Dirichlet  boundary 
conditions  on  some  sides  of  the  polygon  and  Neumann  boundary  conditions  on  others. 
The  results  of  our  pending  investigation  will  be  reported  at  a  later  date. 
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n 

t 

E 

K 

Ti 

180 

0.041426 

9.88098491916389  x  10"15 

4.7357 

r2 

228 

0.086169 

4.13558076672871  x  10"15 

6.1593 

r3 

219 

0.033944 

9.99200722162641  x  10"16 

7.5586 

r4 

289 

0.064797 

3.60822483003176  x  10”15 

18.363 

r5 

704 

1.428900 

5.55111512312578  x  10'16 

21.365 

r6 

1007 

4.284600 

1.22124532708767  x  10”15 

27.666 

Table  1:  Numerical  results  for  the  exterior  Neumann  problem 


n 

t 

E 

K 

r  i 

162 

0.078327 

4.88498130835069  x  10~15 

28.640 

r2 

250 

0.075677 

1.33226762955019  x  10~14 

5.4339 

r3 

303 

0.046949 

3.05311331771918  x  10“16 

31.302 

r4 

349 

0.068338 

5.21804821573824  x  10”15 

107.14 

r5 

776 

0.425080 

2.86437540353290  x  10"14 

99.830 

r6 

997 

0.781490 

1.08801856413265  x  10“14 

46.737 

Table  2:  Numerical  results  for  the  interior  Dirichlet  problem 
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Figure  7:  A  right  triangle  T4  in  R2 


Figure  8:  A  star-shaped  curve  T5  in  R2 
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40 


8  Appendix  A 

In  this  section  we  provide  a  proof  of  Theorem  3.1,  which  is  restated  here  as  Theorem  8.4. 
The  following  lemma  provides  the  value  of  a  certain  contour  integral. 


Lemma  8.1  Suppose  that  —l<p<l  and  0  <  a  <  1  are  real  numbers.  Then 


/  - 7 — \ - o  dz  =  —7 r — t — r 

1  —  2 z  cos(7TQ;)  +  z2  sin(7ra) 


(150) 


Proof.  Let  e  and  R  be  real  numbers  such  that  0  <  e  <  1  and  1  <  R.  Suppose  that  C 


Ce 

-O- 

-e  0  e 


Figure  10:  A  contour  in  C 

is  a  contour  consisting  of  the  intervals  [— R,  —  e]  and  [e,  R],  together  with  the  arcs  Cr, 
defined  by  the  formula  Rel9  for  all  0  <  9  <  n,  and  Ce  defined  by  the  formula  eel6  for  all 
0  <  9  <  7T  (see  figure  10).  We  observe  that 


1  —  2z  cos(7ra)  +  z2  (et1Ta  —  z)(e~ 


(151) 


is  a  meromorphic  function  of  z  inside  C,  with  a  simple  pole  inside  C  at  ema .  Since  e® 
is  within  the  contour,  by  the  residue  theorem 


Jc  1  —  2z  cos(7tq:)  +  z2  Jc  ( eina  —  z){e  mOL  —  i 
We  observe  that  when  for  any  z  G  Cr,  \z\  =  R  and 


dz  =  2m- 


(einay 


ZU 

1  —  2z  cos(7TQ;)  +  z2 

{ei7ra  -  z)(e~i7ra  -  z) 

z)  ~  (R  —  l)2’ 


sm(7ra) 

(152) 


(153) 
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Consequently, 


lim  /  - 7 - - 

R- s-oo  JCr  1  —  2  z  cos(7ro;)  +  z2 


oh 


<  lim 

R-^oo 


'CR 


1  —  2  £  cos(7tq;)  +  z2 


oh 


oh 


<  lim  [  Rft 

R^oo  JCr  (R—  l)2 

Rfj.+  \ 

=  lim  7T— - — - 

R^-oo  (i?  —  l)2 

=  0. 


Likewise,  we  observe  that  for  any  z  G  Ce,  \z\  =  e  and 


ZU 

ZU 

1  —  2 Z  COS(7TO)  +  Z2 

(ei7ra  -  z)(e~i7ra  -  z) 

< 


(i-0: 


(154) 


(155) 


Consequently, 


lim  /  - 7 - - - ; 

e^O  Jc  1  —  2z  COS(7TQ;)  +  Z‘ 


dz 


<  lim 

e->0 


Ce 


1  —  2z  cos(7ra)  +  z2 


dz 


<  lim 


e^°  Jcc  (1  _  e)' 


dz 


rr/4+l 


=  lim  7 r- - ^ 

£->o  (1  -  e)2 

=  0. 


Therefore, 


lim  lim 


i 


dz  = 


R-> oo  e-40  Jc  1  —  2z  cos(7ra)  +  z2  1  —  2z  cos(7ra)  +  z2 

Combining  formulas  (152)  and  (157),  we  observe  that 


(156) 


dz.  (157) 


!  _OQ  1  —  2 z  cos(7TQ;)  +  z2 


dz  =  -71- 


Dll TOtfl 


sm  to 


(158) 


The  following  corollary  evaluates  the  integral  in  formula  (150)  when  it  is  taken  from 
0  to  oo  instead  of  from  —  oo  to  oo. 


Corollary  8.2  Suppose  that  —  1  <  p,  <  1  and  0  <  a  <  1  are  real  numbers.  Then, 


x M  sin(7ra) 


/0  1  —  2x  cos(7ra)  +  x2 


dx  =  n 


sin(/i7r(l  —  a)) 
sin(/i7r) 


(159) 
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Proof.  Lemma  8.1  states  that 


0nrafi 


—  7T- 


sin(7ra)  J_00  1  —  2zcos(7ra)  +  z2 
which  we  rewrite  in  the  form 


—  7 r 


sin(7ra)  J  o  1  —  2zcos(7ra)  +  z1 

poo  poo 

dz  +  ei7Tfl 


dz, 


dz  T 


f_oo  1  —  2z  cos(7ro)  +  Z2 


(160) 


dz 


J  o  1  —  2z  cos(7ra)  +  z2 

Replacing  a  with  1  —  a,  we  also  observe  that 


o  1  +  2z  cos(7TQ;)  +  z2 


—  7T 


gZ7r(l— a)fj,  roo 


sin(7ra)  J o  1  +  2 z  cos(7ra)  +  z2 


dz  +  e47rM 


*o  1  —  2z  cos(7ra)  +  z2 


dz.  (161) 


dz.  (162) 


Multiplying  formula  (161)  by  — e  17111  and  adding  it  to  formula  (162),  we  observe  that 


Di7r(a— l)/i  _  pin(l—a)fi 


7T- 


=  (em^  -  e~mii) 


sin  7ra 


/0  1  —  2 z  cos(7TQ;)  +  z2 


Therefore, 


7T- 


sin(/i7r(l  —  a)) 


sin(7ra) 


sin(/i7r) 


/0  1  —  2z  cos(7TQ;)  +  z2 


dz. 


dz.  (163) 


(164) 


A  simple  analytic  continuation  argument  shows  that  identity  (159)  in  corollary  8.2 
is  also  true  for  all  real  0  <  a  <  2.  This  observation  is  summarized  by  the  following 
theorem. 


Theorem  8.3  Suppose  that  —  1  <  p,  <  1  and  0  <  a  <  2  are  real  numbers.  Then, 


aAsin(7ra)  sin(/i7r(l  —  ct)) 

-  “  —  ClCC  ——  7T  -  ~ 

2x  cos(7TQ;)  +  x 1  sin(//7r) 


(165) 


Proof.  We  observe  that  the  right  and  left  hand  sides  of  identity  (159)  are  both  analytic 
functions  of  a ,  for  all  complex  a  such  that  0  <  Re(a)  <  2.  Therefore,  by  analytic  con¬ 
tinuation  (Theorem  2.9),  it  follows  that  identity  (159)  holds  for  all  complex  a  such  that 
0  <  Re(a)  <2.  ■ 


A  simple  change  of  variables  gives  us  the  following  theorem. 
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Theorem  8.4  Suppose  that  —  1  <  n  <  1  and  0  <  a  <  2  are  real  numbers.  Then 


a^sin^a)  ,  ,  sm(/x7r(l  -  a)) 

- 7 — r - -  ax  =  n  or  - - — 7 — r — - 

2 ax  cos(7ro;)  +  x2  sm(/x7r) 


for  all  a  >  0. 


(166) 


Proof.  Let  a  >  0.  We  observe  that 

a^sin(7ra;) 


J o  1  —  2x  cos(7tq;)  +  x2 

■  „  x^sin (na 

=  a  ^  1  - - — 


dx  = 


sin(7ro;)  dx 


J0  a2  —  2 ax  cos(7tq;)  +  x2 
Therefore,  by  theorem  8.3, 


0  1-2(5)  cos (TO)  +  (|)' 

dx. 


(167) 


x M  sin(7ro;) 

2ax  cos(7tq;)  +  x2 


dx  =  7T 


_j  sin(/i7r(l  —  a)) 
sin(/i7r) 


(168) 
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